
f = figure;                
colormap(f, bone);         
fontname(f, 'Arial');      
ax = axes('Parent', f);    
[X,Y] = meshgrid(alphavec, betavec);

contourf(ax, X, Y, rmat, 20, 'LineStyle','none');
caxis(ax, [0, 2]);
hold(ax, 'on')
plot(ax, alphax, betax, 'p', ...
    'MarkerFaceColor','white', 'MarkerEdgeColor','yellow', 'MarkerSize',13);
hold(ax, 'off')

xlabel(ax, 'Agglomeration in productivity (alpha)', 'FontSize',14);
ylabel(ax, 'Agglomeration in amenities (beta)', 'FontSize',14);

ax.FontSize = 14;
axis(ax, 'square')

c = colorbar(ax);
c.Label.String   = 'Relative Population of Location 1 to 2';
c.Label.FontSize = 14;
